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Abstract. We consider the Renyi entropies S„(£) in the one dimensional spin- 
1/2 Heiscnbcrg XX chain in a magnetic field. The case n = 1 corresponds to the 
von Neumann "entanglement" entropy. Using a combination of methods based on 
the generalized Fishcr-Hartwig conjecture and a recurrence relation connected to 
the Painleve VI differential equation we obtain the asymptotic behaviour, accurate 
to order 0(£~ 3 ), of the Renyi entropies S n (£) for large block lengths £. For 
n = 1,2, 3, 10 this constitutes the 3, 6, 10, 48 leading terms respectively. The o(l) 
contributions arc found to exhibit a rich structure of oscillatory behaviour, which 
we analyze in some detail both for finite n and in the limit n — > oo. 
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1. Introduction 

Let |\&) be the ground state of an extended quantum mechanical system and p = |*) (^1 
its density matrix. In order to quantify the bipartite entanglement in the ground state 
one divides the Hilbcrt space into a part A and its complement B and considers the 
reduced density matrix = Trg p of subsystem A. A measure of the quantum 
entanglement in the ground-state is provided by the Renyi entropies [1] 

S n = ^—lnTrp^. (1) 

The particular case n = 1 of (1) is known as the von Neumann entropy Si and it 
is usually called simply entanglement entropy. However, the knowledge of S n for 
different n characterizes the full spectrum of non-zero eigenvalues of p^ (see e.g. [2]) 
and provides significantly more information on the entanglement than the more widely 
studied von Neumann entropy. 

Of particular interest is the universal scaling behaviour exhibited by S n at 
quantum critical points. For a one-dimensional critical system whose scaling limit 
is described by a conformal field theory (CFT) of central charge c and A being an 
interval of length £ embedded in an infinite system, the asymptotic large-£ behaviour 
of the Renyi entropies is given by [3, 4, 5] 

S n {£)^- Q (l + ^j\nl + c' n . (2) 

Here c' n is a non-universal constant. The scaling behaviour (2) has been verified 
both analytically and numerically for a variety of quantum spin chains whose scaling 
limits are described by CFTs, see e.g. [6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17] 
as well as in direct field theory calculations [18]. In one dimensional systems these 
entanglement entropies provide a very useful way for determining the central charge c 
that characterizes the behaviour at conformally invariant critical points. While other 
methods for determining c such as the finite-size scaling of the ground state energy 
[19, 20] require the knowledge of certain non-universal properties such as the velocity 
of sound, the large-£ behaviour of the entanglement provides a direct measure of c 
as is apparent from Eqn (2). For this reason a scaling analysis of S n is increasingly 
used in numerical studies of quantum phase transitions in one dimensional systems 
[21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33]. 

In such applications S n (£) is computed numerically and the large-^ behaviour 
is then fitted to the form (2). It has been observed that the asymptotic result is 
sometimes obscured by large, and often oscillatory, corrections to scaling [34, 35, 24, 
25]. In Ref. [35], on the basis of both exact and numerical results, it has been 
argued that these corrections are in fact universal and encode information about the 
underlying CFT beyond what is captured by the central charge alone. More precisely, 
they give access to the scaling dimensions of some of the most relevant operators in 
the underlying CFT. This conjecture of Ref. [35] has been recently confirmed by using 
perturbed CFT arguments [36]. 

A precise characterization of the subleading terms in S n (£) is then desirable for 
two reasons. First, the knowledge of their structure will be helpful when using (2) 
to extract the central charge from numerical computations of S n (£). Second, the 
subleading terms can be used to infer the scaling dimensions of certain operators 
in the CFT characterising the quantum critical point. This motivates the present 
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study, in which we significantly extend our recent calculation [35] of the sublcading 
corrections to the Renyi entropies in the XX chain. 

1.1. Spin- 1/2 XX Chain 

The Hamiltonian of the XX model on an infinite one dimensional chain is 

H = - E 2^^i + i + ^t, (3) 

I — — OO 

where af'' y,z are the Pauli matrices at site /. The Jordan- Wigner transformation 

\m<l J 

maps this model to a quadratic Hamiltonian of spinless fermions 

(5) 



/ 1 

H = ~ E c i Q+1 + c '+i Q + 2h ( c ' Q ~ 2 

Z=-oo ^ 



Here /i represents the chemical potential for the spinless fermions c;, which satisfy 
canonical anti-commutation relations {q, cj^} = <5;. m . The Hamiltonian (5) is diagonal 
in momentum space and for \h\ < 1 the ground-state is a partially filled Fermi sea 
with Fermi-momentum 

k F = arccos \h\. (6) 

In the following we will always assume that \h\ < 1 so that we are dealing with a 
gapless theory. 

1.2. Entanglement Entropy of the XX chain: Jin-Korepin Result 

A key result regarding on entanglement measures in the XX chain is due to Jin and 
Korepin [8], who obtained the leading large-£ behaviour of S n . Their result takes the 
form 

S J n K {l) = \{l + ln(2£| sinfc F |) + E n , (7) 

where the constant E n has the integral representation 
r0 ° dt 



E n =[l + 

nl In t 



-) r 

nj Jo 



1/1 | \ 1 e -2r 

' (8) 



1 — n 2 \nsinht/n sinht/ sinht 6 

The objective of our work is to determine the sublcading corrections to S„ K (£) for 
large, finite block lengths I. It is therefore convenient to introduce quantities d n (£) 

d n {£) = S n {£)-S J n K {£), (9) 

to which we will refer throughout. 

The remainder of this paper is organized as follows. For the sake of clarity 
we first present a summary of our results in section 2. When then turn to the 
details of our derivations. In section 3 we briefly review one of our key tools, the 
generalized Fisher- Hartwig conjecture. The latter is used in section 4 to determine all 
"harmonic" corrections to the Renyi entanglement entropies. In order to go beyond 
the generalized Fisher Hartwig conjecture we utilize recent developments related to 
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Random Matrix Theory. These are introduced in section 5 and used to determine 
"non-harmonic" terms in the asymptotic expansion for the von Neumann and Renyi 
entropies in sections 6 and 7 respectively. Comparisons between our analytic expansion 
and numerical results are presented in section 8. 

2. Summary of Results 

This section contains a summary of our results. 

2.1. Renyi Entropies of the XX chain: General Result 
Our full result for d n {£) can be cast in the form 



p,q=l 



2p(2q-l) 



(Qn,q) 



cos(2k F £p) A q sm(2k F p£) 



P 



Lk 



+ 



[5<$e 2 « + h.c] 



T 2 

L k 



where 



1 n + 1 

Zf 285n 3 



Lk 
A q 

Qn,q 



B (n) 



;i5(3n 2 - 7) + (49 - n 2 ) sin 2 k F ) +o(l^ , 



= 2£\smk F \ 



1 + 3 



2q - 1 



cos k F . 



1 V2 ' 2n > 
V(l - 2 Q~ 1 \ 

2q-l 

6n 



(10) 

(11) 
(12) 

(13) 



(5 + 7^)^-15(^ + 1) 



(l + 3^"#)cos(M 



i 2 



(14) 



The leading contribution to d n {£) has already been announced in Ref. [35] and is 
given by 

= 2 COs(2fc F f) kA) -2,n Q + Q[i - mintW]) 

1 — n 



2.2. Renyi Entropies of the XX chain: explicit results for S2{£) and Ss(£) 
In the special cases n — 2 and n — 3 our results read 

Q 2 21 cos{4k F £) 

7Q 2 ^cosk F . Inl ^ i 5 + 3sin 2 fc F 



2Q 2>1 cos(2k F £) 1 
d2{£) = "Li + Zf 



sin(2/c F f) 



64 



+ o(L,T 3 ), (16) 
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and 

, Q 3 ,icos(2k F £) Qi ;1 cos(4fc F £) 4Q 3 .icosfc F . 
d 3 W = + rT^Ts qr 5/3 sm(2A^ 

Q| ! cos(6fc F 15 + 2sin 2 fc F 4 cos(fc F )Q§ j sin(4fc F .) 

+ 243L2 + 3 
Q|. lC os(8fc^) 2g 3 , 1 (lll-62sin 2 (fc F ))cos(2fc F £) , / r _ 3 \ 

+ IZp +— sTIf + ° V M" (17) 

2.5. Renyi Entropies of the XX chain: limit n — > 00 

In the limit n — > 00 infinitely many terms in (10) combine to generate a logarithmic 
contribution, whose general expression is is given in Eq. (70). It assumes a particularly 
simple form at half- filling fc F = it/ 2 

it 2 f 2 e odd , 

^ 241^1-1 .even, (18) 
where 6 = exp(-*(l/2)) w 7.12429. 



2.4. von Neumann Entropy of the XX chain 

In the special case n = 1 corresponding to the von Neumann entropy all oscillating 
contributions to (10) vanish. This explains why it is easier to determine the central 
charge from S\ than from Renyi entropies with n > 2 (this is no longer true in the 
presence of boundaries [34], where it is found that oscillations persist in the limit 
n — > 1). Specializing Eq. (10) to n — 1 we obtain 

Si ~i In* + 4-^(1+ cot A£) . (19) 

In this expression, the £~ 2 power-law behaviour is universal [35, 13]. 



3. Entanglement entropy in the XX model 

Let us return to the spin-1/2 XX model on an infinitely long chain (3). The reduced 
density matrix of a block of I contiguous sites can be expressed as 



(20) 



(21) 



PA = detC exp ]T [Inter 1 -!)]. 
\j,ieA 

where the correlation matrix C has matrix elements 

I sin kp(m — n) 

7r(m — n) 

As a real symmeetric matrix C can be diagonalized by a unitary transformation 

UClfl = 6i m (l + v m )/2. (22) 

This implies that the reduced density matrix pe is uncorrclated in the transformed 
basis. The Renyi entropies can be expressed in terms of the eigenvalues vi as 

i 



s n(£) = ^2 e n {v{) , with e n (x) 



i=i 



1 



In 



1 



+ 



1 



(23) 
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More details about this procedure can be found in e.g. Refs. [6, 16, 37]. We note 
that the above construction refers to the block entanglement of fermionic degrees 
of freedom. However, in the case considered here, the non-locality induced by the 
Jordan- Wigner transformation does not affect the reduced density matrix. In fact, 
it can be seen to mix only spins inside the block. This ceases to be the case when 
two or more disjoint intervals are considered [38, 39] and other techniques need to be 
employed [40] in order to recover CFT predictions [41, 42, 38]. 

The representation (23) is particularly convenient for numerical computations: 
the eigenvalues v m of the I x I correlation matrix C are determined by standard linear 
algebra methods and S n (£) is then computed using Eq. (23). In order to obtain the 
universal behaviour in the limit of large block lengths /-S'cowc follow Ref. [8] . We 
introduce the determinant 

D e (X) = dct ((A + 1)1 - 2C) = dct(G) . (24) 

In the eigenbasis of C the determinant is simply a polynomial of degree I in A with 
zeroes = !,...,£}, i.e. 

i 

D t {\) = H{X-uj). (25) 

This implies that the Rcnyi entropies have the integral representation 

Sn(l) = ^JdXe n (X) d]n ^, (26) 

where the contour of integration encircles the segment [—1,1]. The matrix G is a I x I 
Tocplitz matrix, i.e. its matrix elements depend only on the difference between row 
and column indices 

Gjk = gj-k ■ (27) 

In the theory of Toeplitz matrices an important role is played by the Fourier transform 
9(0) of gi 

m = fj f e" g(0). (28) 

The function g{6) is called symbol and in our case takes the form 

fA + 1 6e[k F ,2n-k F ] 
9{ ' \X-\ 9e [0,k F ]U [27r-fc F ,27r] . 1 ' 

On the interval [0, 27r] the function g(9) has two discontinuities at Q\ = k F and 
8 2 = 2tt — k F . 

3.1. The generalized Fisher- Hartwig conjecture 

The Fisher-Hartwig conjecture [43] gives the asymptotic behaviour of the determinant 
of a Toeplitz matrix in the limit where the dimension I of the matrix becomes large. 
This has been used by Jin and Korepin [8] to derive the leading large I asymptotic 
behaviour of the Renyi entanglement entropies. As stressed in Ref. [8], for the Toeplitz 
matrices defined by the symbol (29) the Fisher-Hartwig conjecture has been proven 
by Basor [44] . 
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In order to employ the Fisher-Hartwig conjecture one needs to express the symbol 
g(9) of a Toeplitz matrix in the form 



9{0) = f(0) Y[ e *M^— sgn(e-0 r )] (2 - 2cos(# - 9 r )) a " 



(30) 



r=l 



where R is an integer, a r , b r and 9 r are constants and f(9) is a smooth function with 
winding number zero. The Fishcr-Hartwig conjecture then states that the large-^ 
asymptotic behaviour of the corresponding Toeplitz determinant is given by 



U =1 



E. 



(31) 



where F[f(9)} = exp(^ J Q W d6 In f (6)) and E is a known function of f{9), a r , b r , 
and 9 r . In our case it is straightforward to express the symbol in the canonical form 
(30). As g(9) has two discontinuities in [0, 2tt) we have R = 2. It is useful to define a 
function 



/?A = ^ln 
2tti 



A + l 
A-l 



where the branch cut of the logarithm is chosen such that 

"A + l" 



— 7r < arg 



A-l 



< 7T. 



= /oe : 



2ibokf 



Inserting the ansatz 

h = -b 2 , ai, 2 = , /(0) = /o = const 

into (30) gives 

' 1 9 e [k F ,2TT - k F ] 

^ e -2nib 2 9 e [0, k F ] U [2tt - k F ,2ir} . 
Comparing (35) to (29) we conclude that we require 

b 2 = (3\ + m , 

where m is an arbitrary integer number. We further identify 

f = (A + l) e - 2ib2kF = (A + l) e -2ifcFm e -2ife F /3 A> 

The integer m labels the different inequivalent representations of the symbol g{8), see 
[44]. In their work Jin and Korepin employed the Fisher-Hartwig conjecture for the 
m = representation and obtained the following result for the large-£ asymptotics of 
Dt (A) [8] 



(32) 

(33) 
(34) 

(35) 

(36) 
(37) 



D J t "(\) 



(A + l) 



A + l 
A-l 



L- 2fj2{x) G 2 (l + h)G 2 (l-p x ). 



(38) 



where = 2^|sinA:F| has been introduced in (11). Inserting (38) into (26) and 
carrying out the integral leads to the result for the asymptotic behaviour of the Renyi 
entropy reported in Eq. (7). Expression (7) provides the leading behaviour of S n (£) for 
large block lengths I. It is the purpose of our work to determine (universal) sublcading 
contributions to (7). This is achieved by noting that for the case when the symbol 
g(9) has several inequivalent representations labelled by an integer m the asymptotics 



Universal corrections to scaling for block entanglement in spin-^ XX chains 8 

of the corresponding Toeplitz determinant is given by the so-called generalized Fisher- 
Hartwig conjecture (gFHC) [44], which reads 

Dt{\) rv^^'e-XKM^'EW. (39) 

m 

In our case, the various parameters in (39) are given by 

l ( m) = ln(/ (ro) ) = ln(A + 1) - 2ik F f3 x - 2ik F m, (40) 
b ^= _ 6 (™)=^ +m , (41) 

where m are integers and 

E (m) = [2 - 2 cos(2fc i ,)]" (,3A+m)2 [G(l + p x + m)G(l - Px ~ m)] 2 . (42) 

Here G(z) is the Barnes G-function [45]. We note that the gFHC has been used to 
determine the large-distance asymptotics of various two-point correlation functions in 
Refs [46, 47]. Important properties of the gFHC in our case are 

(i) The exponential increase is representation independent and governed by the 
exponent 



Re(^ m) ) = Re [ln(A + 1)] - ^ Re 

(ii) The leading oscillatory behaviour depends on the representation and is given by 

Im(4 m) ) = Im [ln(A + 1)] - ^ Im In - 2k F m. (44) 

(iii) The power law correction depends on the representation and is characterized by 
the exponents 

a m = (b[ m) ) 2 + (b^) 2 = -2(l3x + mf. (45) 
The real parts of these exponents are 

Re(a m ) = -2Re(^) - 2m(m + 2Re(/3 A )). (46) 
In conjunction with the inequality —1 < 2Re((3\)) < 1 this establishes that 

Re(a m ) < Re(a ). (47) 

Equality in (47) holds only for m = 1 and Re (/3a) = — \, which corresponds to 
the case — 1 < A < 1. 

We note that point (iii) is crucial for Eqn (38) to give the correct asymptotic 
behaviour of S n (£): along the integration contour in (26) we always have Im(A) 7^ 0. 
Representations with m/0 therefore give subleading corrections, which we are going 
to analyze in the following section. 

The full result of the generalized Fisher-Hartwig conjecture for the Toeplitz 
determinant takes the form 

De ^ (X + l yf^±±\ * J2 L -2(m+0 x f e -2 ikF me 

x [G(m + 1 + px)G(l - m - /3 A )] 2 . (48) 



In 


"A + 


1" 










A- 


1 





(43) 
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4. Corrections to the scaling for entanglement 



The leading corrections to scaling for the Renyi entropies are obtained from the 
"harmonic" terms given by the generalized Fisher-Hartwig conjecture. It follows from 
(46) that the most important corrections arise from the first two contributions with 
m = ±1. Keeping only the three terms corresponding to m = —1,0,1 in (48) we 
obtain the following expression for the asymptotics of the determinant D((X) 

_ 2 _ 4fev G 2 (2 + /3 A )G 2 (-/? A ) 



JK 



1 



-2ik F £ 



L 



+e 2ik F e L - 



k 

2+4/3 x 



G 2 (1 + /3 A )G 2 (1-/3 A ) 
G 2 (2-/? A )G 2 (/? A ) 
G 2 (1 + /3 A )G 2 (1-/? A ) 



(49) 



Here D^ K (X) is given in Eq. (38). Using G(l + x)/G{x) = T(x) we can rewrite the 
last formula as 



*/(A) 



= c -2ik F l L -2(l+2l3 x ) T 2 (l + f3\) c 2»fc F lj-2(l-2ftQ r 2 (l - /3 A ) 



(50) 



r 2 (-/? A ) 1 " ^ k r 2 (/? A ) 

It follows from the factorized form of (50) that the contributions of the correction 
terms to the entropies are easier to calculate than the contribution of the leading term 
Dj K itself. This will enable us to obtain a full analytic answer. For large Lk we have 
(we recall that d n (£) = S n (£) - S-{ K (£)) 



d n (£) 



1 

2jd 



dX e„(A) 



d!n[l + */(A)] 



' ' dX 



(51) 



1 

dX 2iri 

The contour integral can be written as the sum of two contributions infinitesimally 
above and below the interval [—1,1] respectively, i.e. 



d n (£) 



1 

2iti 



/1+ie r-l — ie 

-1+ie J — 1—ie 



dX e„(A) 



dX ■ 



(52) 



This shows that we only require the discontinuity across the branch cut. The only 
discontinuous function is j3\, which for —1 < x < 1 behaves as 



Px±ie = ~iw(x) T 



with 



1 1 + x 

w(x) — — In 

w 2tt 1 - x 



We now change variables from A to w 

— oo < w < oo. 

We have 



t-2-4,0 r 2 (i 



A = tanh(7ru;) 



(53) 
(54) 



2 _ 4/3 r 2 (i + /3) 
fe r 2 (-/3) 



^4ito„,2 



r -2+4,3 



r 2 (i-/3) 
r 2 (/?) 



r -2+4,3 



r 2 (i-/3) 
r 2 (/3) 



^j=-iu>+i 



, 7 M, 

7 2 (-w), 



r— 4iu> 2/ 



where we have dropped terms of order 0(L k 4 ) compared to the leading ones and we 
have defined 



T(l -iw) 

nh+iwY 



(55) 
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Integrating by parts and using 
d 



, e„(tanh(7ru;)) = — 
aw 1 



-(tanh(ri7rw) — tanh(7rw)) . 



(56) 



we arrive at 



d n (£) 



in 



2(1 -n) 



f 

j — i 



dw(tanh(nw) — tanh(n7rw)) x 



x [e^^LfVM - e 2ikFe L^ 4iw ^ 2 (-w)] + . . . 



(57) 



For large I the leading contribution to the integral arises from the poles closest to the 
real axis. These are located at w = i/2n (w = —i/2n) for the first (second) term in 
(57). Evaluating their contributions to the integral gives 

2 



1 — n 



+ o(r 2 ' n ). 



(58) 



L I 2 ' In > 
^(2 _ 2n)- 

This result implies that at half-filling (k F — |) and n > 1 the corrections are positive 
(negative) for odd (even) I. 



4.1. Subleading Corrections 

Eqn (58) describes the asymptotic behaviour in the limit L/, — > 00, n fixed. It provides 
a good approximation for large, finite I as long as ln(Lfe) 3> n. This is a strong 
restriction already for moderate values of n. For example, is required to be larger 
than 10 4 for n = 10. For practical purposes it is useful to know the corrections to 
S n (£) for large I but ln(Lfc) not necessarily much larger than n. In this regime there 
are two main sources of corrections to (58). 

(i) The integral (56) is no longer dominated by the poles closest to the real axis 
and contributions from further poles need to be included. These give rise to 
corrections proportional to L k 2q ^ n , with q integer. 

(ii) Further terms in the expansion of the logarithm in Eqn (51) need to be taken 
into account. The corresponding contributions are proportional to e ±l2pfcF ^ with 
p = 2,3,.... 

At half-filling (zero magnetic field) the situation is different in that terms with 
odd p all give rise to an overall factor (— l) e and hence modify the staggered 
contribution to S n (£), while terms with even p add to the smooth (non-oscillatory) 
part already present in S^ K (£). 

We now take both types of corrections into account. We first consider the series 
expansion of the logarithm in Eq. (51) 

h[1+ „ w]= f: '-')' +1 (w . (69) 

p=i 

Recalling the explicit expression (50) for \E^(A) leads to a binomial sum 

= E ( P )e 2ifc ^ (29 - p) ife 2p ife 4( ^ 29)/3A ^7cl^ , (60) 
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where we have introduced the shorthand notation cp = (T(l + 0)/T(— f3)) 2 . When 
calculating the discontinuity across the branch cut running from A = —1 to A = 1 all 
terms other than q = and q = p give rise to terms that are subleading in Lk ■ Hence 
we may approximate 

(^(tanh(™) + ie)) p - (^(tanhfTru;) - ie)) p « e-™°* t *I% v "<?_ iw _ 1/i 

_|_ e 2ik F £p j^-iiwp^j} 



-iw+1/2 ■ 



(61) 



The analog of (57) then reads 



~ (_1)p+i 
P =i 



p 2(1 -n) 



dw(tanh(7rw) — tanh(wra;)) 



(62) 



The integral is carried out by contour integration, taking the two terms in square 
brackets into account separately. The first (second) contribution has simple poles in 



the upper (lower) half plane at w q 
integer such that 2q — 1 ^ n, 3n, 5n, 

d n (£) = 2 



- i 2 -± 



(w q 



-§=-), where q is a positive 



2n — "In 

Contour integration then gives 



1 



~ (-DP+ 1 2 p (2 g -l) 

^ ^ cos{2k F lp)L k » (Q n ,,) 



l-2/n 



), 



(63) 



p,g=l 

where the constans Q„ jg have been defined in (13). In the sum over q the special values 
2q — 1 7^ n, 3n, 5n, . . . are to be omitted. In particular, this means that for n = 1 all 
these corrections are absent. Eqn (63) is one of the main results of our work. It shows 
that there are contributions to the Renyi entropies with oscillation frequencies that 
are arbitrary multiples of 2k F . 

At half- filling (k F = tt/2) certain simplifications occur. For even I we find 



d n (£) 



1 



while for odd i 

d n (l) ~ — — 
1 — n 



(2l)"Q n ,i - (2*)" 



we obtain 



(2*)~ Q„,i + (2*)- 



?n,l 



(2£)" 



+ (2^)" 



<9r 



Qn.l 



Qn,3 



+ Q 



n,3 



+ 



(64) 



(65) 



In all the above analysis we have ignored contributions to the generalized Fisher- 
Hartwig conjecture with \m\ > 1. While these lead to oscillatory contributions with 
frequencies that are integer multiples of 2k F they are suppressed by additional powers 
of i~ x and hence are subeading, even in the case where n is not small. 

It is apparent from (63) that the limit n — > oo deserves special attention. Soo(£) is 
known in the literature as single copy entanglement [48] . Here it is neccessary to sum 
up an infinite number of contributions in order to extract the large-^ asymptotics. We 
note that the large- n limit is not only of academic interest, but will provide information 
on the behaviour of S n (£) in the regime n 3> hiL^, 3> 1. 



4-2. Large n limit of S n (£) 

In order to investigate the limit n — > oo we consider cqn (62), but now first take the 
parameter n to infinity and then carry out the resulting integrals. This gives 

(-l)f 



OO 



P =i 



P 



I 

J — ( 



dw(sgn(w) — tanh(7™)) 
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X [ e -Kk F ep L 4iwp ^ w )fP e 2ik F ep L -Hwp [ 7 („ W )]2P 

00 ( _-i \p r poo 
= -J2 [ -^L e - 2 ^Im / dw[l - tanh(irw)]L 4 k lwp [ 7 H] 2 



P =i 



l2p 



„2ik F £ 



poo 

p lm / dw[l - t&nh(irw)]L- 4twp [j(-w) 
Jo 



2p 



(66) 



Using that the first singularity in the upper (lower) half plane occurs at w = i/2 
(w = —i/2) we deform the contours to run parallel to the real axis with imaginary 
parts i/A and — i/4 respectively, i.e. for the first term we use 



poo /'«/4 poo+i/4 

/ dw f(w) — dw f(w) + dw f(w). 

JO Jo Ji/4 



It is straightforward to show that the second integral contributes only to order 0{1/L k ) 
and does not give rise to logarthimic corrections. Hence the leading contribution is of 
the form 

,1/4 



(-1)" 



P =i 



P 



-2ik, 



ip Rt 



• / dz(l 
Jo 



it<m(irz))L- 4zp ( 7 (w)) 



,2p 



+e 2 lkF i PRc f dz{ i + i t an (7rz))Lfc 4zp ( 7 (iz)) 
Jo 

= 2 V ^ '- cos(2k F £p) / dze- 4zplnLk { 1 {iz)f p . 



(67) 



For large Lk the dominant contribution to integral is obtained by expanding ( 7 (iz)) 2p 
in a power series around z = 

00 ( _^p ,1/4 

doo(*)~2V^ L cos(2k F £p) / dze' 4zplnLk {l + 4pz^{l/2) + ...) 



= *E 



(-1)p 

cos(2k F £p) 

pi 



1 , *(l/2) 
41nL fe 4 In 2 L fr 



+ 0^), 



(68) 



where ^(x) is the digamma function. The leading contribution can be expressed in 
terms of the dilogarithm function Li 2 (x) using 



*E 

P =i 



(-ir 



cos(2fc F £p) = Li 2 (-e i2fep£ ) + Li 2 (-e-" fc ^). 



-i2k F £\ 



P 



(69) 



In the half-filled case (k F = w/2) our result takes a particularly simple form 

1 ^ 2 

(_l)p(*+i) 



— E 



21nL fc ^-; 

P =i 



= < 



21n£ fc 6 
1 tt 2 



£ odd , 



(70) 



£ even . 



21ni fc 12 

Summing some of the subleading terms in (68) to all orders in (ln(L fc )) _1 leads to an 
expression of the form 



Odd 



241n(6L fe ) -1 £even, 



(71) 
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where b = exp(— \I>(l/2)) w 7.12429. This is found to be in good agreement with 
numerical computations. 

5. Connection with random matrix theory 

Keating and Mezzadri [49] have shown that the Toeplitz determinant D((X) is related 
to an important quantity in random matrix theory, namely the gap probability for 
the circular unitary ensemble (CUE). The generating function Ef VE [(0, 0); £] (in the 
following we will drop all arguments to ease notations) for the probability of finding 
exactly k eigenvalues e l9 within the segment 9 e (tt — (ft, ir] of the unit circle is given 
by [50] 

*< CUE - (£ -* D * ■ ■ ■ (£ -< /' ) *,i - 



This is equal to the determinant of the Toeplitz matrix 



i<j<k<e 



ir_i=0, a; = l, zi = --; 1— . (75) 



= tu W) with Wn = 8 nfi + -^(-1)" +1 . (72) 

£71% ft 

It then follows that for £ = 2/(A + 1) and (ft = 2/c F we have [49] 

D e (\) = (\ + l) e E? lJE . (73) 

For any value of £ the generating function Ef VE can be determined from a recurrence 
relations connected to the Painleve VI transcendent [50]. The recurrence relation 
reads 

xtxt-i - c = [(£ + l)x e+1 + {£- l)x t -i] - 1 ~ Xi ~ 1 [tx t + (£- 2)x e - 2 ] , (74) 

where c = cos kp and the initial values are 

£ sin kp 

Kl-^kp 

The generating function is related to xe by 

Ef+i'Ef-i' _ 1 2 _ Di+iDj-i , 7a -. 
{e cve )2 - L x e~ D 2 ■ W 

For the sake of completeness we quote the values of the generating function for £ = 
and £ = 1 

£ CUE = 1, E? VE = l--^<ft. (77) 

27T 

5.i. Leading large-£ asymptotics of xg 

In Ref. [49] it was suggested to combine the asymptotic results (38) following from the 
Fisher-Hartwig conjecture with the recurrence relation (74) in order to obtain further 
corrections to the large-^ behaviour of the Renyi entropies. Inserting (38) into (76) 
suggests that [49] 

Xe= V2\M +0{g -2y (?8) 
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However, a numerical solution ol the recurrence relation shows that (78) does not 
generally provide the correct large-£ asymptotics of xg. The reason for this is as 
follows. When we substitute the "full" result (50) of the generalized Fisher-Hartwig 
conjecture into (74) we find that the contributions due to the representations with 
m = ±1 behave as £~ 1±2 P>- for large I. For any Re(p\) ^ one of these will dominate 
over the contribution arising from the m = term that gives rise to Eq. (78). In other 
words subleading contributions to Dg{\) give rise to the leading large-£ behaviour of 
x e \ 

We now show in more detail how to extract the large-£ behaviour of xg from that 
of Dg{\). In order to keep things simple, we focus on the case Re(p\) > 0. Here we 
may neglect the terms with \m\ > 1 and m = — I in (48), which leads to 



D 2 

where we have introduced 

a = (2 sin ftp) 



Mi 



[1 - Aa 2 e 2 ^ e r 2+4f3 >~ sin 2 k F ) + . . . , (79) 



-1+2/3. 



(80) 



The contribution in square brackets arises from the m — Fisher-Hartwig term and is 
the result quoted in Rcf. [49]. For Re(p\) > this term is subleading and we obtain 
instead 



x t ~ (-l)V feF ¥ 



{2smk F ) 



Re(p x ) > 0. (81) 



Here we have fixed the sign of xg by requiring that the expression (81) asymptotically 
satisfies the recurrence relation (74). The analogous analysis in the case Re((3\) < 
gives 

T(1 + /3a) 



xg - (-l) l e' lkFl i- 1 - 2h {2sink F y 



2[3x 



Re(/? A )<0. (82) 



,..(83) 



r(-&) 

We may combine Eqns (81) and (82) into a single equation 

x t ~t¥ [e*^sinM^ r T ^+ e ^^^ 

« L 1 (Pa) 1 (-Pa) 

We emphasize that (83) must not be understood as giving the two leading terms in 
the large-£ asymptotic expansion of xg because e.g. for Re(p\) > g there are other 
contributions to xg that decay more slowly than i^ 1 ^ 2!5x . In order to check the result 
(83) we have solved the recurrence relation (74) numerically for a number of different 
values of A and k F - In Fig. 1 we compare the asymptotic expression (83) against the 
numerically computed values for xg. The agreement is seen to be excellent in all cases. 



5.2. Asymptotic expansion for xg and analytic corrections to the gFHC expression 
forDg(X) 

We now turn to the derivation of contributions to the large-£ asymptotic expansion 
for Dg{\) that are not contained in the gFHC. This will be achieved by utlizing the 
recurrence relation (74). 

We expect the asymptotics of Dg{\) to be such that each harmonic term predicted 
by gFHC is multiplied by a function analytic in 1/L Restricting out attention to the 
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Figure 1. Real part of xg as a function of I for several values of A and kp. 
The points are obtained from a numerical solution of the recurrence relation 
(74). The continuous lines are obtained from the asymptotic prediction (83) 
by the replacement (— 1) — > e 17rt . The first two panels correspond to the same 
value p x 2t 0.323792 - 0.128075i but two different values of k F . The last panel 
corresponds to P\ ~ 0.0158924 — 0.0003966i and hence the contributions of both 
terms in Eq. (83) are important. In all cases we observe good agreement of the 
theoretical prediction (83) with the numerical data. 



first three harmonic terms (i.e. m — —1, 0, 1), this leads to an expansion of the form 
Pi 



1 



Cl 

I 



C2 

P 



_j_ a 2 e 2ik F eg-2(l-2P x ) 



1 



at 



«2 
P 



+ b 2 e -2ik F i^-2(l+2p x ) 

where do is given by (80) and 



6l bn 

1 + — + — 
£ P 



(84) 



b a = (2sin&p) 



-l-2/3 A 



r(i + /?A) 



(85) 
(86) 



We note that by definition we have 

(2 sin k F ) 2 a Q b Q = . 

We now proceed in a straightforward albeit extremely tedious way: 

(i) We first insert (84) into (76) in order to obtain an expression for the asymptotic 
expansion for x^. 
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(ii) Wc input the resulting expression into the recurrence relation (74) for xg and 
determine the parameters characterizing the asymptotic expansion of X£ order by 
order in £~ . 



The result of step (i) for Re (/3a) > is 

3=0 



asy 



1 Ojl Oj2 Oj3 



a 2j e 2ik F j£ p(-2+4p x ) 



+ 



3=1 



£-4j/3x e -2ik F je 

(2ao sin /cf) 2 -? 



i 



qj2 

e 2 



+ .... 



(87) 



where we have written only the terms required for our purposes and where have 
introduced the quantity 

xf = (-l)V fc ^r 1+2 ^2a sinfc F . (88) 

The explicit expressions for the coefficients Oji and qji in terms of the expansion 
coefficients aj , bj , and Cj characterizing the large-f asymptotics of Dp are reported in 
Appendix A. 

Step (ii) consists of substituting (87) in the recurrence relation (74) and 
determining the coefficients Oji and q 3 \. The non-linearity of (74) renders this a very 
difficult task, because terms at different orders in I" 1 in xp, contribute to the same 
order in the recurrence relation. For this reason it is crucial to retain sufficiently many 
terms in (87). The results of this procedure are reported in Appendix B. 

Combining the results reported in Appendix A and Appendix B then yields the 
desired expressions for the expansion coefficients aj , bj and Cj 

Cl (/3 A ) =20licatk F , 

c 2 (/3a) = y ("I + 7Pl + 12 Pt - 3/3a( 5 + ^l) csc " k ?) > 



6,(/3 A ) = Cj (/3 x + l) , j = l,2, 
a j {0x) = c j {0 x -l), j = 1,2. 



(89) 



6. Corrections to the von Neumann entropy 

Having determined the asymptotic expansion for De(\) we may now use (26) to 
calculate additional subleading contributions to the Renyi entropies. We first consider 
the von Neumann entropy (the case n = 1), in which as we have seen above all 
harmonic contributions vanish. Taking the limit n — > 1 in (26) and following through 
the same steps as in section 4 we find 



di(£) 



f 

J — c 



dw- 



TTW 



cosh TTW 
where we have defined 



c, 1 — c 



l _|_ 2c t 



(c+) 2 -2c 2 - + (cr) 2 
2P 



4 



Cj(—iw =F 



1, 



1,2. 



Here Ci j2 are given by (89) and we have used 

tanh(7rw) — tanh(n7rw) ttw 

1 — n cosh 2 ttw 



lim 



(90) 



(91) 



(92) 
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As c+ — c 1 is an even function of w the 0(1 x ) contribution in (90) vanishes. The 
0(£~ 2 ) contribution can be calculate analytically using the integrals 



dw- 



irw 



COsh TTW 



f 

J — ( 



dw- 



nw 



cosh TTW 



7 
120 



which gives the final result 

di(*)~- 



12^ 2 



+ cot kp 



(93) 



(94) 



The simplicity of this answer suggests the existence of a much more straightforward 
derivation than ours. 



7. Corrections to the Renyi entropies 



The case of the Renyi entropies is more complicated because the contribution of the 
harmonic terms does not vanish. Our starting point is the expansion (84) for the 
Tocplitz determinant, which we express in the form 



Pi 

D JK 



(95) 



Here ^>i{j3\) arc the contributions we have taken into account previously in section 4. 
The logarithm of the Tocplitz determinant is expanded as 



hi 



r D t ' 




~ In 






OO 




P =i 



1 + . f(ft)+ »2Ws) + 



M\/3 x )[Hf t (l3 x )] 



p-i 



(!)/ 



1 2 



1P-2 



X (0) (/?a) + 



X (1) (/?a) , X (2) (/?a) 



(96) 



Following through the same steps as in section 4 we then arrive at the following 
expansion for the Renyi entropies 

coo 



d n (£) 



in 



2(1 -n) 



dw (tanh(7rw) — tanh(n7ru>)) 



x (o) (ft0 + 



X (1) (/?a) , X i2 HM 



+ 



£ 2 



&\=-iw-\ 



= d^{l) + d^{t) + d^{l). (97) 

The contribution dffl {i) has been determined in section 4. The other two contributions 
are calculated by the same method as in section 4 and we find 



2 cos(/c_f 
1 - n 



J2 (-l) p+1 sin(2fc F p£) L] 



1 2p(2g-l) 



p,q=l 
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2p 



1 + 3 



2q-l 



2 _ 2 _2p(M-V L 

— J2(-l) P L k 2 - 

p,q=l 



pi 1 i 2 q-l 
1 2 T 2n 



r /l _ 2g-l 
-I [ \ 2 2n 

r(§ + ^) 
V 2 2n J . 



18 



(98) 



2p 



B(») e 2<PfcF< + hx- 



l^^z + l/ ^ 15(3^-7) 

+ ^ 1440n 3 ^ 71 + gin 2 fcF 

Explicit expressions for the coefficients are given in Appendix C. 



(99) 



8. Numerical results 



Given our asymptotic expansion a natural question to ask is how well it approximates 
the Renyi entropies for large but finite block lengths I. In order to address this 
question we will now present a number of comparisons between our asymptotic result 
and numerically exact expressions for S n (£). The latter are obtained by determining 
the eigenvalues of the Toeplitz matrix C nm in Eq. (21) and computing S n from Eq. 
(23). 



8.1. Leading contributions to d n (£) 

In the top two panels of Fig. 2 we plot the absolute value of d n (£) for n — 2 and 
n = oo at kp = it/2 and compare it to the leading asymptotic expressions (15) and 
(71) respectively. We see that the asymptotic expressions give good agreement with 
the numerically exact results even for moderate values of £. 

The next issue we turn to is the behaviour of S n (£) for large, finite values of n. 
In this case the asymptotic power-law d n (t) cx l^ 2 l n only emerges for very large block 
lengths ln^ ^> n. On the other hand, for smaller values of I the numerical data is 
seen to follow the large n prediction (71) as is shown in the bottom panel of Fig. 2. 
Here we plot l/\d n {t)\, which at large values of I will grow as l 2 l n . A logarithmic 
behaviour for small small values of i is clearly visible, which then crosses over to the 
expected £ 2 ^ n regime at approximately ln^ ~ n. It appears that the crossover scale 
is larger for even £. We expect a crossover between these two regimes to be a generic 
feature in critical theories. This suggests that in such gapless models particular care 
is required when studying S n (£) for large n. 



8.2. Corrections to the von Neumann entropy S\{£) 

For the von Neumann entropy S\ all the oscillating terms vanish and the predicted 
large-£ asymptotic behaviour is given in Eq. (19). In Fig. 3 we plot —£ 2 d\{£) as a 
function of £ and compare it with the prediction (19). We see that the agreement 
is excellent, which indicates that further corrections are very small. We note that 
for vanishing magnetic field (fcp = 7r/2) the amplitude of the 0(£~ 2 ) correction term 
is numerically small (1/60) so that the corresponding contribution to S\{£) becomes 
negligible already for relatively small £. At least in the particular case of the XX 
model in zero field this shows that the central charge is most conveniently extracted 
from finite-size scaling studies of Si (£) rather than higher Renyi entropies. 
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/ 

Figure 2. Top: d n (l) = S n (£) - S^ K (£) at half-filling (k F = f ) for n = 2 and 
n = oo compared to the the asymptotic expressions (straight lines for even/odd £ 
respectively). The agreement is seen to be excellent even for moderate values of 
I. Bottom: |d n (f)| _1 as a function of I for several values of n and kp = 2, The 
straight lines show the asymptotic results (71) in the limit n — > oo for even and odd 
I respectively. We see that for large n the correction d n {£) exhibits a logarithmic 
increase up to a block size ln£ ~ n, when the asymptotic behaviour starts to be 
seen (as we are plotting |d n (^)| _1 the asymptotic behaviour corresponds to a £ 2 / n 
power- law increase with (). 

8.3. Subleading contributions to d n (i) 

For n > 1 the structure subleading corrections to scaling for S n (£) is significantly 
richer. Explicit expressions, accurate to order 0(£~ 3 ), for the cases n — 2 and n = 3 
are given in Eqns (16) and (17) respectively. A comparison of these results for n = 2 
to numerical computations is presented in Fig. 4. We have chosen kp = 7r/4 so as to be 
able to separate the oscillation frequencies of the various contributions. The top panel 
in Fig. 4 presents a comparison of the asymptotic expression for d,2(£) (continuous 
lines) to numerical computations (dots) and show good agreement even for small £. In 
order to better assess the accuracy of the asymptotic results we introduce the rescaled, 
subtracted quantity 

D 2 {£) = [d 2 {£)-d^{l)]£\ (100) 

where d^ y {£) represents the leading correction given in Eq. (15). By construction 
Z?2 (£) should tend to a sum of oscillatory terms with fixed amplitudes for large £. The 
four-sublattice oscillatory behaviour of D^it) predicted by (16) is clearly visible and 
as expected we observe excellent agreement between the numerical and asymptotic 
results. For larger values of n, the corrections arising from the 'analytic' part of Di(X) 
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Figure 3. Correction to scaling for the von Neumann entropy di = S± — S^ K . 
We plot the quantity — £ 2 di(£) for values kp = nir/20, where n = 1 (top curve) 
up to n = 10 (bottom curve). The straight lines are our prediction Eq. (19). 
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Figure 4. Corrections to scaling for the Renyi entropy S2(£) at kp = 7r/4. 
Upper panel: d2 (t) as a function of i. Dots are the numerical results while 
the continuous line corresponds to the asymptotic expression (16). Lower panel: 
rescaled subleading corrections D2(£) defined in (100) as a function of i. The 
agreement between the asymptotic expression (continuous line) and numerical 
data (dots) confirms that (16) is correct to order o(£~ 2 ). 



are less important than the harmonic contributions. Thus the most relevant terms 
in the asymptotic expansion are those given in Eq. (63). In fact, the first analytic 
correction has an exponent — 1— 2/n and at a given n appears only after [n/2] harmonic 
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contributions with exponents —2p/n with p = 1, 2, . . .. It has already been observed 
in Ref. [35] that for higher values of n the first subleading order does not suffice do 
give an accurate description of the Renyi entropies. 

In Fig. 5 we show a comparison of the corrections d n (£) for n — 10,20 and 
kp = 7r/4 with the asymptotic result Eq. (63). Step by step we take into account 
further terms in the asymptotic expression (63) until we obtain good agreement with 
the numerical data. We observe that for n — 10 three terms in Eq. (63) are enough to 
reproduces the data, while for n = 20 we need five terms to have the same accuracy. 




/ 



Figure 5. Corrections d n (£) for n = 10 (top) and n = 20 (bottom) for kp = it/4. 
The numerical data are well described by Eq. (63), but more terms are needed to 
obtain the same degree of accuracy when n is increased. In both plots the various 
continuous curves correspond to Eq. (63) one (red curve), two (green curve), 
three (blue curve), etc terms in (63) retained. 



9. Conclusions 

In this work we have determined the asymptotic behaviour of the Renyi entropies S n (£) 
in the spin-1 /2 XX model for large block lengths I. A summary of our results has been 
presented in section 2. While we have considered the specific case of the spin-1/2 XX 
chain in a magnetic field, some features we find are in fact universal. In particular, the 
scaling of the leading oscillatory term (15) has been observed for the XXZ model in zero 
magnetic field in Ref. [35]. The corresponding exponent is modified to £- 2K / n j where 
K the Luttinger liquid parameter. This is in full agreement with recent perturbed 
CFT calculations [36] . As we have emphasized repeatedly, a precise knowledge of the 
structure of the oscillating terms in S n {^£) is useful for extracting properties such as the 
central charge and scaling dimensions of certain operators at quantum critical points. 
They furthermore can be used for analyzing numerical studies of more complicated 
quantities such as the entanglement of two disjoint intervals [38, 40]. 

Oscillating behaviour has also been observed in numerical studies of other 
entanglement estimators [51] such as the valence-bond entanglement. A natural 



Universal corrections to scaling for block entanglement in spin-^ XX chains 22 

question is whether these can be determined for certain models using the free fermion 
techniques we employed for the XX case as well. 

Finally we would like to remark that our results carry over directly to the critical 
Ising chain. According to Ref. [12], the Renyi entropies in the critical Ising model 
(with c = 1/2) are related to those of the spin-1/2 XX chain in zero magnetic field 
{kp = tt/2) by 

S%(e) = ±SZ x (2£,k F = */2). (101) 
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Appendix A. Relation between the expansion coefficients for Di and xi 



In this appendix we report the relations between the coefficients (Zj, 6j and c, in 
the asymptotic expansion (84) of the Toeplitz determinant Di(\) and the expansion 
coefficients Oji, qji characterizing the large-£ behaviour (87) of the auxiliary quantities 
X£. The following relations hold: 



Oji = i(l — 2j3\) cot k F + 



oi - Ci 



(2j + 1) , 



Oji = U + g )(«2 - c 2 ) + (1 - 3f3 x + 3pl) + i cot k F (2j + 1 - (j + l)0)(ai - ci) 



+ 2j ^ L (ai - ci)((2j - l)ai - (2j + 3) Cl ) - fy, 



2 sin kp 



910 = - 

911 = Pi 



3. 

Oi - Ci 



+ - 2/3 A )cot k F 



_ (02 - c 2 )/3^ _ ^_ ^ a lCl (^ + 2) 3a?$ ^ c?($ + 8) _ $(1 - 2/3 A + 5$) 



9i2 = „ - 3c 2 + 



+ ((oi - Cl )$l + (1 - 20 A )ci)t cot fc F + 



$(3-14& + 15$) 



4 sin fcp 



920 = . 

921 = £ A 



ai + &i 



4^cotfc F -^(^ + l) Cl . 



(A.l) 



We have derived relations for some other coefficients such as 003 and 013, but they are 
not needed for our purposes and we therefore refrain from reporting them here. 



Appendix B. Expansion coefficients for 

Substituting the expansion (87) into the recurrence relation (74) gives a set of 
consistency relations for the coefficients qji and Oji in (87). These read 

920 = 930 = 940 = 921 = 931 = 941 = , (B.l) 
on = 3o i - 2(1 - 2/3 x )i cot kp , (B.2) 
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9n 

Of)2 



gio(-ooi +2/3 x icotk F ) , 
ooi + 2.7(001 - (1 - 2/3 A )« cot k F ) 
/3 A (-l + 12/3 A -32/3| + 27^) 
6 

/3 A (2-13/3 A 



+ 



32^ 
4 sin 2 kp 



9 12=^O02-^(4 



9^) 



1M) 
2 sin 2 ftp 



(13 + 18/3 2 ) . 



(B.3) 
(B.4) 



(B.5) 
(B.6) 



Appendix C. Expressions for the coefficients B pq 

We recall the expressions for the coefficients ci,2 and 61,2 (89) 
Cl {l3x) = 2pli cot k F , 

c 2 (ftO = ^(-1 + 7/3^ + 12/31 _ 3; g2 (5 + 4/3 2 } csc 2 fcF) 



6,(/3 A ) = Cj -(l + &) , j = 1,2. 
,1 

2g- 1 1 



In terms of the constants £>+ , c+ g 



2n 



C i,Q = C i 



2q-l 1 
2n 2 



(C.l) 

J = 1,2 • (C.2) 



the coefficients B pq are given by 



B$ = 2sin 2 (fc F ) 



2,q 2 2 ><? 

If 2q — 1 = n, in, 5n, ... we instead have Bp™g = 0. An explicit expression is 
B$ = * [(5 + 28x 2 ) sin 2 (M - 15(4x 2 + 1)] 

- P -[(l + l2x 2 )cos(k F )\ 



(C.3) 



(C4) 
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